Mon. Not. R. Astron. Soc. OOP. [TlfT6l C2002') Printed 6 June 2011 (MN WT^ style file v2.2) 



Thermal emission at 4.5 and 8 [im of WASP-17b, an 
extremely large planet in a slightly eccentric orbit 

D. R. Anderson/* A. M. S. Smith/ A. A. Lanotte/ T. S. Barman,^ 
O A. Collier Cameron/ C. J. Campo/ M. Gillon/ J. Harrington/ C. Hellier/ 
<^ P. F. L. Maxted/ D. Queloz/ A. H. M. J. Triaud/ P. J. Wheatley^ 

^ t ^Astrophysics Group, Keele University, Staffordshire ST5 5BG, UK 

^ ' ^ Institut d'Astrophysique et de Geophysique, Universite de Liege, AUee du 6 Aout, 11, Bat. BSC, Liege 1, Belgium 
■ ^Lowell Observatory, I4OO West Mars Hill Road, Flagstaff, AZ 86001, USA 

^SUPA, School of Physics and Astronomy, University of St. Andrews, North Haugh, Fife, KY16 9SS, UK 
^Planetary Sciences Group, Department of Physics, University of Central Florida, Orlando, FL 32816-2385, USA 
^ Observatoire de Geneve, Universite de Geneve, 51 Chemin des Maillettes, 1290 Sauverny, Switzerland 
"^Department of Physics, University of Warwick, Coventry, CV4 7AL, UK 

Accepted 1988 December 15. Received 1988 December 14; in original form 1988 October 11 

ABSTRACT 

We report the detection of thermal emission at 4.5 and 8 |J.m from the planet WASP- 
17b. We used Spitzer to measure the system brightness at each wavelength during two 
occultations of the planet by its host star. By combining the resulting light curves 
with existing transit light curves and radial velocity measurements in a simultaneous 
analysis, we find the radius of WASP-17b to be 2.0 i?jup, which is 0.2 i?jup larger than 
any other known planet and 0.7 -Rjup larger than predicted by the standard cooling 
theory of irradiated gas giant planets. We find the retrograde orbit of WASP-17b to 
be slightly eccentric, with 0.0012 < e < 0.070 (3 a). Such a low eccentricity suggests 
that, under current models, tidal heating alone could not have bloated the planet 
to its current size, so the radius of WASP-17b is currently unexplained. From the 
measured planet-star flux-density ratios we infer 4.5 and 8 |^m brightness temperatures 
of 1881 ± 50 K and 1580 ± 150 K, respectively, consistent with a low-albedo planet 
that efficiently redistributes heat from its day side to its night side. 

Key words: methods: data analysis - techniques: photometric - occultations - plan- 
ets and satellites: atmospheres - planets and satellites: individual: WASP-17b - stars: 
individual: WASP-17 ~ planetary systems - infrared: planetary systems. 



1 INTRODUCTION 

WASP-17b is a transiting, 0.49-Jupiter-mass planet in a 
3.74-day orbit around a metal-poor, 1.3-Solar-mass F6V star 
l|Anderson et al.ll2010bl. hereafte r AlO) . Data presented in 
AlO, confirrn e d by iTriaud et all (|2010l . hereafter TIO) and 
iBayliss et all l|2010i '). indicated that WASP- 17b is in a ret- 
rograde orbit around its host star, the first such orbit to be 
found. 

WASP- 17b probably formed beyond the ice line in a 
near-circular, coplanar orbit, as predicted by the canon- 
ical model of star and planet formation. It may have 
subsequently acquired an eccentric, highly inclined or- 



bit via planet-planet scatteri ng (e.g. [ Ford &: Rasiol 



or the Kozai mechanism (e.g. iFabryckv fc Tremaind 



2008) 



2007 
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iNagasawa. Ida fc Besshj I2OO8I ). Tidal friction may then 
have shortened and circularised the long, eccentric orbit, 
with the energy being dissipated with in the planet as heat 
(e.g., Jackson. Greenberg fc Barneai2008bl : llbgui fc Burrow^ 
120091 1 

As the planet is low-mass and the star is hot (Tcff — 
6650 ± 80K, TIO) and fast-rotating (usin/ = 9.8 kms~\ 
TIO), the radial velocity (RV) measurements are relatively 
low signal-to-noise, and so the eccentricity of the planet's 
orbit was poorly constrained (e — 0-0.24) in the discovery 
paper (AlO). This translated into uncertainties in the stellar 
radius {R*— 1.2-1.6 Rq) and thus in the planetary radius 
(-Rpi — 1.5-2.0 -Rjup), meaning the planet is larger than pre- 
dicted by the standard cooling theory of irradiated gas giant 
planets by 0.2-0.7 R.^up (Fortney, Marlov fc Barnes 2007). 

Using a coupled radius-orbit evolutionary model, 
llbgui fc Burrowsl (|2009l ) demonstrated that planets can be 
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inflated to radii of 2 -Rjup and beyond during a transient 
phase of heating that accompanies the tidal circularisation 
of a short (a « 0.1 AU), highly eccentric (e ~ 0.8) orbit. 
Such a large radius persists for only a few tens of Myr, sug- 
gesting that we are unlikely to observe any one system dur- 
ing this brief stage. However, only 3-4 of the ~100 known 
planets are extremely bloated and the transit technique does 
pref erentially flnd large planets. 

iLeconte et al.l (,2010) argued that the tidal heating rate 
is underestimated for ev en moderately ec c entric orbits in 
studies such as that of llbgui fc Burrowi (|2009l ) . If true, 
then a large fraction of energy tidally dissipated within the 
planet would have been radiated away by the age typical 
of the most bloated planets (a few Gyr) and so could not 
have played a significant role in their observed bloating. 
llbgui. Sp iegel fc Burrow"^ (|201ll ) admit that their equations 
are not applicable at large eccentric i ty, bu t counter that nei- 
ther are those that iLeconte et al.l (|2010l ) use. They state 
that the current uncertainty in tidal theory means that no 
approach can be considered correct. 

In AlO, the derived radius of WASP-17b is largest 
(2.0 i?jup) when a circular orbit is imposed, smaller 
(1.7 -Rjup) when eccentricity is a free parameter (e = 
0.13), and smaller again (1.5 -Rjup) when a main- 
sequence prior is imposed on the star an d ecce ntric- 
ity is let free (e = 0.24). llbgui fc Burrows! (|2009h and 
llbgui, Spiegel fc Burro^ (|201ll ) each note that, compared 
to planets that did not undergo tidal heating, tidally in- 
flated planets are still signicantly larger Gyr after their or- 
bits have circularised and tidal heating has ceased. In each 
study though, the orbits are still significantly non-circular 
(e > 0.1) when the planets are largest. Hence, if the or- 
bit of WASP-17b were found to be near-circular (which 
would mean that Rpi « 2 iijup) then it would seem unlikely 
that the planet could have been inflated by a single episode 
of tidal inflation. Another possi bility is an ongoing tidal 
heatin g scenario, as explored by llbgui. Burrows fc Spiegell 
l|2010h . in which the orbit of WASP-17b would be kept non- 
circular by the continuing interaction with an as-yet undis- 
covered third body. 

In order to better constrain the stellar and planetary 
radii, the system age, and the potential transient and ongo- 
ing tidal heating rates, an improved determination of orbital 
eccentricity is required. Further high-precision RV measure- 
ments were obtained and presented in TIO. These allowed 
eccentricity to be better constrained to e = 0.0661q'q43- 
best prospect of improving this situation further lay with the 
measurement of an occultation (as the planet passes behind 
the star), which would constrain ecosw, where u is the ar- 
gument of periastron. 

In addi tion, by observing o ccultations from the 
ground (e.g. 'Anderso n et al.l l2010al ) and from space (e.g. 
[Whcatlcv ct al. 2010), we are able to perform photometry 
and emission spectroscopy of exoplanets which are spatially 
unresolved from their host stars. This allows us to deter- 
mine planet albedos and the rates at which energy is redis- 
tributed from the day side to the night side of the planet (e.g. 
IBarman 200§,), and to infer the temperature structure (e.g. 
Knutson et al.ll2008[) and chemica l composition of planet at- 



mospheres (e.g. Swain et al.ll2009l '). 

We present here observations of two occultations of 
WASP-17b, each of which was measured at both 4.5 and 



8 \im. We combine these new data with existing data in a 
simultaneous analysis to show that WASP-17b is the largest 
{Rpi — 2.0 -Rjup) and least-dense (ppi = 0.06 pjup) planet 
known, and is in a slightly eccentric orbit around a 2-3 
Gyr-old, F-type star. Exoplanet occultation photometry is 
at the limit of Spitzer systematics and reliable conclusions 
concerning atmospheres and orbits depend on a careful anal- 
ysis of the data. We thus present a detailed description of 
our method. 



2 NEW OBSERVATIONS 

We observed two occultations of the planet WASP-17b by 
its Ka = 10.22 hos t star with the Spitzer Space Telescope 
(|Werner et al.1 1200-1') during 2009 April 24 and 2009 May 1. 
On each date we measured the WASP-17 system with the 
Infrared Array Camera (IRAC. [Fazio et aLll2004l l in full ar- 
ray mode (256 x 256 pixels, 1.2 arcsec/pixel) simultaneously 
in channel 2 (4.5 (xm) and channel 4 (8 ^.m) for a duration 
of 8.4 hr. We used an effective integration time of 10.4 s, 
resulting in 2319 images per channel per occultation. To 
redu ce the time-depende nt sensitivity of the 8 [im channel 
(e.g. iKnutson et al.ll2008l l. we exposed the array to a bright, 
diffuse source (M42) for 214 frames immediately prior to 
each occultation observation. 

We used the images calibrated by the standard Spitzer 
pipeline (version S18.7.0) and delivered to the community as 
Basic Calibrated Data (BCD). We added back to each image 
the estimated brig htness of the zodiacal light in the sky dark 
(jReach et al.1120051 ) so photometric uncertainties and the op- 
timal aperture radii are correctly determined. For each im- 
age, we converted flux from MJy sr~^ to electrons and then 
used IRAF to perform aperture photometry for WASP-17, 
using circular apertures with a range of radii (1-6 pixels). 
The apertures were centred by fltting a Gaussian profile on 
the target. The sky background was measured in an annu- 
lus extending from 8 to 12 pixels from the aperture centre, 
and was subtracted from the flux measured within the on- 
source apertures. We estimated the photometric uncertainty 
as the quadrature addition of the uncertainty in the sky 
background in the on-source aperture, the read-out noise, 
and the Poisson noise of the total background-subtracted 
counts within the on-source aperture. We calculated the 
mid-exposure times in the HJD (UTC) time system from 
the MHJD_OBS header values, which are the start times 
of the DCEs (Data Collection Events), by adding half the 
integration time (FRAMTIME) values. 

We found (see Section 3) that for WASP-17 the highest 
signal-to-noise is obtained when using an aperture radius of 
2.9 pixels for the 4.5 |i.m data, and a radius of 1.6 pixels for 
the 8 [im data. The data are displayed raw and binned in 
the flrst and second panels, respectively, of Figure [1] 

We rejected any flux measurement that was discrepant 
with the median of its 20 neighbours (a window width of 
4.4 min) by more than four times its theoretical error bar. 
We also performed a rejection on target position. For each 
image and for the x and y detector coordinates separately, 
we computed the difference between the fitted target posi- 
tion and the median of its 20 neighbours. For each dataset, 
we then calculated the standard deviation a of these me- 
dian differences and rejected any points discrepant by more 
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Table 1. Number of points rejected per dataset per criterion. 



Dataset 


Flux 


x-pos 


y-pos 


Total (per cent) 


2009 Apr 24 / 4.5 


34 


4 


1 


35 (1.5) 


2009 Apr 24 / 8 n-m 


29 


31 


30 


45 (1.9) 


2009 May 1 / 4.5 ^m 


35 


6 


6 


40 (1.7) 


2009 May 1 / 8 H-m 


17 


25 


22 


37 (1.6) 



than 4 a. The numbers of points rejected on flux and target 
position for each dataset are given in Table [1] 

According to the IRAC handbook, each IRAC array 
receives approximately 1.5 solar-proton and cosmic-ray hits 
per second, with ~2 pixels affected in channel 2, and ~6 
pixels per hit affected in channel 4, while the cosmic ray flux 
varies randomly by up to a factor of a few over time-scales 
of minutes. Thus the probability per exposure that pixels 
within the stellar aperture will be affected by a cosmic ray 
hit is 1.5 per cent for channel 2 and 1.3 per cent for channel 
4, which is in good agreement with the small portion of 
frames that we rejected. 



3 DATA ANALYSIS 
3.1 Data and model 

We determined the system parameters from a simultaneous 
analysis incorporating: our new Spitzer occultation photom- 
etry; the WASP discovery photometry covering the full or- 
bit for the three seasons (March to August) of 2006-2008 
and presented in AlO; a high-precision, /c-band transit light 
curve taken with the 1.2-m Euler-Swiss telescope on 2008 
May 6 and presented in AlO; and 124 RV measurements, 
including 34 taken during transit, made with the CORALIE 
and HARPS spectrographs and presented in AlO and TIO. 

These data were input into an adaptive Markov-chain 
Monte Carlo (MC MC) algorithm ('Colli er Cameron et al] 
l2007l : |Pollacco et al. , 2008 ; Enoch et al. 20101). Such a simul- 
taneous analysis is necessary to take account of the cross- 
dependency of system parameters and to make an hon- 
est assessment of their uncertainties. We used the follow- 
ing as MCMC proposal parameters: Tc, P, A_F, Tn, b, K\, 
Teff, [Fe/H], y^coscj, y/esinoj, \/ v sin / cos A, V v sin J sin A, 
A^4.5(i.m and A_fg^m (see Table [5] for definitions). 

At each step in the MCMC procedure, each proposal 
parameter is perturbed from its previous value by a small, 
random amount. Stellar density, which is constrained by the 
shap e of the transit light curve (jSeager fc Mallen-Ornelaa 
|2003| ) and the eccentricity of the orbit, is calculated from 
the proposal parameter values. This is input, together with 
the latest values of Tofi and [Fe/H] (which are controlled 
by Gaussian priors ) into the empirical mass calibration of 
lEnoch et all l|2010l ) to obtain an estimate of the stellar mass 
Mf From the proposal parameters, model light and RV 
curves are generated and is calculated from their com- 
parison with the data. A step is accepted if (our merit 
function) is lower than for the previous step, and a step 
with higher x^ is accepted with probability exp(— Ax^/2). 
In this way, the parameter space around the optimum so- 
lution is thoroughly explored. The value and uncertainty of 



each parameter are respectively taken as the median and 
central 68.3 per cent confidence interval of the parameter's 
marginalised posterior probability distribution. 

As Ford ( 2 006 ) notes, it is convenient to use ecosu and 
esinoj as MCMC proposal parameters, because these two 
quantities are nearly orthogonal and their joint probabil- 
ity density function is well-behaved when the eccentricity 
is small and uj is highly uncertain. Ford cautions, however, 
that the use of ecoscj and esinoj as proposal parameters 
implicitly imposes a prior on the eccentricity that increases 
linearly with e. As such, we instead use y^costj and y^sinoj 
as proposal parameters, which restores a uniform prior on e. 
For similar reasons, we use V v sin / cos A and V v sin I sin A 
rather than v sin I cos A and v si n I sin A to parameteri se the 
Rossiter-McLaughlin effect (e.g. ICaudi fc WinnllioOTi ). 



3.2 Spitzer data 

3.2.1 Deciding between models and datasets 

Systematics are present in IRAC photometry at a level simi- 
lar to the predicted planetary occultation signal. Therefore, 
it is necessary to carefully detrend the photometry so as 
to obtain accurate occultation depths and timings. To dis- 
criminate between various detrendin g models we ti sed the 
Bayesian Information Criterion (BIC. ISchwar j[l978l ) : 



BIC = x +fclnA^ 



(1) 



where k is the number of free model parameters and A'^ is 
the number of data-points. The BIC prefers simpler models 
unless the addition of extra terms significantly improves the 
fit. As such, it is a useful tool for selecting between models 
with different numbers of free parameters. 

In addition, we used the root mean square (RMS) of the 
residuals about the best-fitting trend and occultation models 
to discriminate between different light curves obtained from 
different reductions of the same sets of images. 



3.2.2 Aperture radii 

We determined the optimal aperture radii to use for the 4.5 
and 8 \im data by performing aperture photometry with a 
range of aperture radii (1-6 pixels), and choosing the radii 
that produced the maximal signal-to-noise (Figure [2l e.g. 
.Howell 1989 ). 

For the 4.5 [im data, the radius that results in the high- 
est signal-to-noise is 2.9 pixels and it is this radius that we 
adopt (Figure [21 upper panel). This radius incorporates the 
majority (~92.5 per cent) of the target flux and little back- 
ground flux (~2.6 per cent of that of the target). 

At 8 |i.m, as compared to 4.5 (xm, the background is 
brighter by a factor ~20 and the source is fainter by a fac- 
tor ~4. Thus, with increasing aperture radius, the back- 
ground flux quickly dominates the target flux (Figure [2J 
lower panel). Indeed, the background flux equals the tar- 
get flux when using our adopted, optimal (highest signal- 
to-noise) radius of only 1.6 pixels. This radius incorporates 
~60 per cent of the target flux and a similar amount (97 per 
cent of the target) of background flux. The background flux 
is greater than the target flux by factors of approximately 
1.5, 2, 3 and 7.5 within apertures with radii of 2.2, 2.75, 3.5 
and 6 pixels respectively. 
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Figure 1. In each of the above three plots, the upper two datasets were obtained at 4.5 |xm (blue) and 8 ]im (red) on 2009 April 24 
and the lower two datasets were taken at 4.5 |J.m (blue) and 8 ]im (red) on 2009 May 1. Relative flux offsets were applied to datasets 
for clarity. Left: Raw Spitzer data with the best-fitting trend and occultation models superimposed. Middle: The same data binned in 
phase (A0 = 0.002, ~11 min) with the best-fitting trend models superimposed. Right: The binned data after dividing by the best-fitting 
trend models, and with the best-fitting occultation models superimposed. We normalise the flux received from the star alone to unity, 
which is measured during occultation. 



3.2.3 Systematics 

IRAC uses an InSb detector to detect light around 4.5 vim, 
and the measured flux exhibits a strong correlation with the 
position of the target star on the array. This efi^ect is due to 
the inhomogeneous intra-pixel sensitivity of the detector and 
is well-documented (e.g.iKnutson et al.ii200^. and references 
therein) . Following ICharbonneau et al.l ( 20081 ) we modelled 
this effect as a quadratic function of the sub-pixel position 
of the PSF centre, but with the addition of a linear term in 
time: 



Table 2. Trend model parameters and coefficients 



df = ao + a^dx + Uydy + 



xdx^ + Uyydy'^ + atdt 



(2) 



where df = f ~ f is the stellar flux relative to its weighted 
mean, dx — x — x and dy = y — y are the coordinates 
of the PSF centre relative to their weighted means, dt is 
the time since the beginning of the observation, and ao, 
Qx, ay, a^x, CLyy and at are coefficients. We determined the 
trend model coefficients by linear least-squares minimiza- 
tion at each MCMC step, subsequent to division of the data 
by t he eclipse rnodel. We used singular value decomposi- 
tion l|Press et al.lll993 ) for this purpose. Though a common 
eclipse model was fitted to occultation data from the same 
channel, trend models were fitted separately to each dataset. 

The best-fitting trend models are superimposed on the 
binned photometry in the middle panel (first and third 
curves from the top) of Figure [1] Table [2] gives the best- 



4.5 i^m 

2009 Apr 24 2009 May 1 



8 [un 

2009 Apr 24 2009 May 1 



ao 



106955.09 
24.28 
25.28 
-26.0t^j 
-2606.7ti33 
-5049.5^20 3 

207.7l^g;2 



107395.47 
24.85 
25.44 

-io.8ii°:l 

4927 O"*"®^'^ 

-68m.5tlll 

''-eei.g 
3663.4j:}^5.7 



17289.49 
25.23 
23.02 
-108. 42^ 



-105. 0^ 



,+0.67 
-0.77 

720.2^6.2 

42.3l^-^ 

-1368.7111*? 

-136:0 
-134.3 

648.9t|-2 



-1246. S 



17866.66 
25.64 
23.18 
-94.2 
-142.0 
-1301 
-1748.4 
905.4 
557.7t|;° 



+3.1 
3.1 
+ 19.1 
-20.1 
+7.7 
-7.2 
+401.7 
398.4 
122.2 
121.2 



fitting values for the trend model parameters (Equation [2]) , 
together with their 1-a uncertainties. 

We found consistent 4.5 [im eclipse depths when in- 
corporating one of the two datasets or both of them in 
our analysis: AF^.s = 0.00225 ± 0.00015, AF4.5 pun., = 
0.00244 ± 0.00020, and AF4.5 \im,^^ = 0.00230 ± 0.00012. 

The systematics in the data from 2009 April 24 are of 
much smaller amplitude than those in the data from 2009 
May 1. This is due to a chance placement of the target star 
on the detector. The detector positions of WASP-17's PSF 
centre are shown in Figure [S] During each occultation and 
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Figure 2. The flux due to WASP-17 (red circles), the sky and 
instrumental background (green triangles), and both WASP-17 
and the background combined (blue squares) , as well as the signal- 
to-noise (magenta diamonds) as a function of aperture radius. 
The upper and lower panels show the 4.5 |im and 8 |^m data 
respectively. We show the data from 2009 Apr 24, but the data 
from 2009 May 1 produce near-identical plots. 



in each channel the motion due to the nodding of the space- 
craft is evident in the x and y positions of the PSF centres. 
However, contrasting the two occultations, there is a marked 
difference in the radial distance from the nearest pixel centre 
over the course of the observations. In the data from 2009 
Apr 24 the placement of the target on the detectors is such 
that the motion of the spacecraft in the x-direction largely 
compensates for the motion in the y-direction, resulting in a 
near constant radial distance from the nearest pixel centre. 
The opposite is the case in the 2009 May 1 dataset, where 
the motion in the x and y directions combines to produce 
large-amplitude oscillations in the distance from the nearest 
pixel centre. This results in the large saw-tooth systematics 
seen in the hght curve (second panel of Figure [1]). 

Near the beginning of the observations on 2009 May 1, 
the target crossed a pixel boundary on the 4.5 yim detector 
(Figure O middle-right plot). This resulted in a point of 
inflection in the distance of the target from the nearest pixel 
centre (Figure [S] bottom-right plot). As the sensitivity is 
higher toward the pixel centre and lower near the edges, it 
is therefore curious that no corresponding inflection point is 
seen in the hght curve (Figure [1] middle panel). 

IRAC uses a Si As detector to observe at 8 [im, and 
its response is usually thought to be homogeneous, though 
another systematic affects the photometry. This effect is 
known as the 'ramp' because it causes the gain to in- 
crease asymptotically over time for every pixel, with an 
amplitude dependin g on a pixel's illumination history (e.g. 
iKnutson et al.|[2008l . and references therein). Again follow- 



ing [Charbonnea^^et|al] l|2008l ). we modelled this ramp as a 
quadratic function of \n{dt): 

df = ao + ai ln{dt + tos) + a2{\n{dt + toff))^ (3) 

where toff is a proposal parameter. To prevent toff from drift- 
ing more than an hour or so prior to the first observation, we 
place on it a Gaussian prior by adding a Bayesian penalty 
to our merit function (x^)- 



BPt., 



,2 / 2 
toff/o-t„ 



(4) 



where at^^ = 15 min. 

From an initial MCMC run, we observed systematics 
in the residuals of the second 8 (xm dataset, and so investi- 
gated decorrelating the 8 |i.m data with detector position. A 
significantly lower occultation BIG (ABIC — —141) resulted 
when also detrending for detector position, i.e. detrending 
with Equation [2] rather than with Equation |3l In addition, 
there was less scatter in the 8 \im data when decorrelating 
with detector position (i.e. when detrending with Equation^ 
rather than with Equation [S] Figure |4] and Table [Sjl . When 
not decorrelating with detector position, significantly dif- 
ferent best-fitting 8 \im occultation depths were obtained 
for the two individual datasets (Table [Sjl and the depth ob- 
tained from the combined datasets was much deeper than 
otherwise. For these reasons and for reasons that will be pre- 
sented in the remainder of this section, we opted to decor- 
relate the 8 [im data with detector position. 

In Section r3.3l we use deconvolution photometry to show 
that the observed dependence on detector position is likely 
to have been introduced during aperture photometry. There- 
fore there is no evidence of an inhomogeneous intrapixel re- 
sponse of the SiAs 8 [im detector, contrary to the case with 
the InSb 4.5 \im detector. 

The second and fourth curves from the top in the middle 
panel of Figure [1] are the best-fitting trend models when 
detrending the two 8 [im datasets with Equation [2] Note 
that the saw-tooth patterns of the 8 \im trend models are in 
phase with those of the 4.5 [im trend models, though each 
dataset was fit separately with its own trend model. 

In addition to Equation [51 which we will call lin- 
ear TIME, we tried trend functions with a variety of time- 
dependency. These were NO time: 



df — spatial 



(5) 



where spatial = ao + a^dx + aydy + a^xdx^ + ayydy^ repre- 
sents the detector position terms and an offset; QUAD time: 



df — spatial + atdt + attdt 



linear ln time: 



df — spatial + ai \n{dt + toft); 

QUAD LN TIME: 

df = spatial + ai ln(dt -I- toff) -I- a2(ln(cit -I- tos))^; 
and RISING EXP (|Harrington et al.ll2007l ): 
df — spatial + 03 exp(a4dt). 



(6) 



(7) 



(8) 



(9) 



where 04 is a proposal parameter. 

In Table |4] we present the occultation depths and BIG 
values resulting from detrending the 8 |i.m data with the 
various models. The 4.5 and 8 |i.m data were fitted simul- 
taneously, so an improved fit to the 8 yuai data would not 
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Figure 3. The detector positions of WASP-17's PSF centre for the first occultation (left-hand plots) and the second occultation (right- 
hand plots). The PSF centre positions on the 4.5 and 8 |im detectors are respectively depicted by blue and red dots. For each occultation 
we show the distance of the PSF centre from the nearest pixel centre in the x and y directions (top and middle panels respectively) and 
in the radial direction (bottom panels). Pixel centres are located at (x, y) = (0, 0) while pixel edges are located at (0.5, 0.5) and are 
demarcated by dashed lines. 



Table 3. A comparison of the 8 \aa occultation depths and residuals from deconvolution photometry and aperture photometry. 



Method 


Trend cq. 


AF8 ^ml+2 


A-Fg i^mj 


A-Fg |xm2 


RMS8|xmi+2 


RMSg^mi 


RMSg^m^ 


aper. phot. 


El 


0.00238 ± 0.00036 


0.00187 ± 0.00050 


0.00303 ± 0.00059 


0.01103 


0.01142 


0.01010 


aper. phot. 


O 


0.00455 ± 0.00036 


0.00246 ± 0.00058 


0.00643 ± 0.00055 


0.01112 


0.01127 


0.01119 


decon. phot. 


O 


0.00240 ± 0.00036 


0.00202 ± 0.00052 


0.00276 ± 0.00047 


0.00994 


0.01019 


0.00967 



be preferred if the fit to the 4.5 \im data were considerably 
worse. In Table|4]the models are presented in descending or- 
der by how well they fit the combined 8 \im datasets, and the 
BIG values are given relative to the best-fitting model (lin- 
ear time). The LINEAR TIME model is strongly favoured and 
we thus adopt this as our trend model for the 8 \x.m data. 
The BIG values resulting from MGMG runs incorporating 
the two 8 |J.m datasets individually (8vimj and 8vtm2) are 
also given, with a similar order of preference to that of the 
combined datasets. Again, the linear time model is clearly 
favoured over the others, supporting our decision to use the 
same model for the two datasets. The occultation depths are 



consistent between the four most preferred models, but are 
not so for the two less preferred models. The depths found 
from the first dataset are shallower than those found from 
the second dataset, with a difference between the two of 1.7a 
in the case of the linear time trend model. 

In Figure [5] we present the detector positions of WASP- 
17 both during and outside of occultation. On 2009 Apr 24, 
the star occupied the same region of the detector during the 
occultation as when outside of occultation. However, on 2009 
May 1, the star occupied different regions of the detector 
during occultation than when outside of occultation, though 
there was some overlap. The reason for this can be seen 
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Figure 4. A comparison of the methods for reducing and detrending the 8 \^.m data. In each of the above three panels, the upper 
three curves (red dots) are the data from 2009 Apr 24 and the lower three curves (blue dots) are the data from 2009 May 1. In each 
triplet of curves, the top curve is the light curve obtained by aperture photometry and detrended v^ith Equation |3] The middle curve 
is the light curve obtained by aperture photometry and detrended with Equation [2] The bottom curve is the light curve obtained by 
deconvolution photometry and detrended with Equation[3] Left panel: Binned raw data, with the best-fitting trend models superimposed. 
Middle panel: Binned detrended data, with the best-fitting occultation models superimposed. Right panel: Binned residuals about the 
best-fitting trend and occultation models. 



in the top-right panel of Figure [31 which shows that the 
star moved steadily in the x-direction. This was in addition 
to the motion due to the nodding of the spacecraft, which 
resulted in some overlap between the in-occultation and out- 
of-occultation detector positions. As we decorrelate the light 
curves with detector position, the data from 2009 Apr 24, 
with the greater detector position overlap, are thought to be 
more reliable. However, the data from the two occultations 
detrend similarly well, and we find no reason to disregard 
the latter dataset. 



This indicates that, though we had requested the same 
detector positions for the target for each observation run, 
small differences in the pointing and motion of Spitzer can 
result in markedly different systematics. 

We investigated using the fitted PSF positions from the 
higher signal-to-noise 4.5 [im data in the aperture photome- 
try and positional decorrelation of the 8 [im data. To account 
for the offset between the two detectors we fit the differences 
in the x and y directions and translated the coordinates by 
those amounts. The 8 yim occultation depths, both when in- 
corporating one of the two datasets or both of them, were 



8 D. R. Anderson et al. 



Table 4. The 8-|J.m occultation depths and the combined (4.5 and 8 |a.m) relative occultation BIC values, when detrending the 8 nm 
data with the various models 



Model 


Eq. 


AFg |xmi+2 




AFg i^nij 


ABICg^m,+, 


ABICg|j.mj 


ABICgixmj 


linear time^ 


m 


0.00238 ± 0.00036 


0.00187 ± 0.00045 


0.00307 ± 0.00055 











quad In time 


m 


0.00271 ± 0.00042 


0.00226 ± 0.00060 


0.00329 ± 0.00063 


14 


38 


31 


quad time 


m 


0.00291 ± 0.00059 


0.00179 ± 0.00085 


0.00397 ± 0.00086 


16 


17 


15 


rising exp 


m 


0.00260 ± 0.00037 


0.00188 ± 0.00047 


0.00307 ± 0.00055 


17 


17 


17 


linear In time 


m 


0.00363 ± 0.00038 


0.00318 ± 0.00046 


0.00417 ± 0.00054 


60 


34 


25 


no time 


m 


0.00210 ± 0.00037 


0.00103 ± 0.00046 


0.00353 ± 0.00058 


252 


210 


18 



tPor LINEAR TIME: BIC4.5n.mi^2+8Hmi+2 = 10 772, BIC4.5Hmi^2+8|xmi = 8 298 and BlCi.syim^^^+syxm^ = 8 072. 



very similar to those obtained when fitting the stellar PSF 
position in the 8 [im data, and there was no significant re- 
duction in the residual scatter about the best-fitting models. 
Hence, we proceeded as before. 

Spitzer's pointing oscillates around the nominal posi- 
tion, with an amplitude of ~0.1 pixels over a period of ~1 
hour. We also see higher frequency jitter, with periods of 1-2 
minutes (the cadence of our data is 1 2 seconds), in the posi - 
tion of WASP-17. Some authors (e.g. IWheatlev et al.ll201(]| ') 
chose to smooth the measured target positions prior to fight 
curve detrending. However, we found that detrending with 
the unsmoothed positions resulted in a reduced BIC (ABIC 
— —931), and in smaller residual RMS values: 5.3 and 8.8 
per cent lower for the two 4.5 \im datasets, and 0.6 and 1.4 
per cent lower for the two 8 \im datasets. 

To ascertain whether the observed short-period jitter 
was due to measurement error, we measured the position of a 
second star in the field for the two 4.5 (im datasets. For both 
WASP-17 and the second star we subtracted their Gaussian- 
smoothed (cr = 84 s) positions to remove the longer-period 
oscillations. We then fitted Gaussians to the distributions of 
the detector x and y coordinates of both stars and of their 
relative separations. If the measured positions of WASP-17 
and the second star are uncorrelated, then the variance of 
the distribution of relative separations would be the sum of 
the variances of the distributions of each star's positions. 
However, we found that the distribution of separation in 
the i-direction had a variance smaller than that by a factor 
nine for the first dataset and by a factor two for the second 
dataset. For the y-coordinate, the factors were 25 and 6 for 
the two datasets. Thus, the short-period jitter is real and the 
light curves should be detrended with unsmoothed target 
positions. 



3.2.4 Aperture radii revisited 

As a check of the choice of aperture radius (2.9 pixels) for the 
4.5 [im data, we input the 4.5 vim light curves obtained with 
each aperture radius into a simultaneous MCMC analysis 
that incorporated all but the 8 ^im data. These analyses 
produced consistent 4.5 \im occultation depths (Figure [6l 
upper panel), indicating that the 4.5 jtm result is relatively 
insensitive to the choice of aperture radius. 

As a check of the choice of aperture radius (1.6 pix- 
els) for the 8 ^im data, we input each 8 ^tm fight curve 
into a simultaneous MCMC analysis that incorporated all 
other data. When decorrelating with detector position (Fig- 
ure|6l middle panel), the fitted 8 yuai occultation depth varies 



weakly with aperture radius. Beyond an aperture radius of 
3.5 pixels (by which point the fiux due to the sky background 
is 3 times that of the target within the target aperture), a 
deeper occultation is measured. Without decorrelating with 
detector position (Figure |S1 lower panel) , the fitted 8 [im 
occultation depth is a strong function of aperture radius. 

As the conclusions drawn from Spitzer occultation 
observations depend on accurately measured occultation 
depths, we advise others to check for a correlation between 
fiux and detector position in their 5.8 [im and 8 [im datasets, 
and for a dependence of occultation depth on aperture ra- 
dius. For example, from Figure 1 of .Fressin et al.. (,20101 it 
appears that similar patterns of saw-tooth systematics are 
present in both the 4.5 and 8 [im fight curve, though they 
only decorrelate the former fight curve with detector posi- 
tion. If a dependence of measured fiux on detector position 
was introduced during aperture photometry then the mea- 
sured 8 (i.m occultation depth could be erroneous. 



3.3 Deconvolution photometry 

To verify the 8 |i.m occultation depths and to investigate 
the source of the dependence of the measured 8 [im fiux 
on detector position, we obtained 8 vim fight curves by 
performing deconvolution photometr y with d e cpho t . This 
method was first described by iGillon et al.l (|2006l . |2007| ) 
and has been optimized for Spitzer data by Lanotte et al. 
(in prep). It is based on th e imag e-deconvol ution method 
of Magain. Courbin fc Sohvl l| 19981 . see also iMagain et al. 
(2007^), which respects the sampling theorem of IShannonI 
(1949), in contrast with traditional deconvolution methods. 
In a first step, 25 random BCD images taken on 2009 April 
24 were used to determine a partial PSF. This was then used 
to deconvolve the whole set of images and to determine op- 
timally the position and fiux of WASP-17. 

The DECPHOT light curves do not exhibit a position- 
dependent modulation of the fiux (Figure|3)). Therefore, the 
saw-tooth modulation seen in the light curves obtained from 
aperture photometry (Figure |4]) is likely due to a pixellation 
effect rather than an intra-pixel inhomogeneity in IRAC's 8 
[im detector response. During aperture photometry of the 
8 vim data, an aperture radius of only 1.6 pixels was used. 
The calculation of a circular aperture is non-trivial and the 
majority of photometry routines make a polygonal approx- 
imation, which tends to be less accurate for smaller radii. 
Aside from that the calculation of how much flux should be 
attributed to partial pixels is another potential source of er- 
ror. A better result is obtained if a PSF is used, rather than 
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Figure 5. The detector positions of WASP-17 both during and 
outside of occultation. The top panel shows the occultation of 
2009 Apr 24 and the bottom panel shows the occultation of 
2009 May 1. The 4.5 (blue squares) and 8 (red circles) [xm data 
taken outside of occultation are overplotted with the 4.5 (magenta 
saltires) and 8 (cyan crosses) |j.m data taken during occultation. 
The 4.5 ]im data are shown relative to detector position {x,y) 
= (24,25), and the 8 \im relative to {x,y) = (25,23). Note that, 
though the axes' ranges are the same between the two plots, each 
abscissa covers only 60 per cent the range of each ordinate. 



if uniform illumination is assumed, but even that is not per- 
fect. Partial deconvolution is a photometric method that is 
optimal in a least-squares sense, i.e. the background contri- 
bution is minimized because each pixel is properly weighted. 
As this is not the case for aperture photometry, and as the 
background at 8 \im is bright relative to the target, we had 
to use a small aperture to optimize the signal-to-noise ra- 
tio of our measurements, leading to pixelisation effects that 
translated into a correlation of the measured flux with de- 
tector position. 
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Figure 6. Top panel: The dependence on aperture radius of 
the fitted occultation depth (red circles with error bars) and the 
residuals (blue up-triangles = 2009 Apr 24 data, green down- 
triangles = 2009 May data) for the 4.5 \ira data. Middle panel: 
The same as the top panel, but for the 8 ]im data and when 
treating the 'pixel phase' effect. Lower panel: The same as the 
top panel, but for the 8 [ira data and when neglecting the 'pixel 
phase' effect. 



We performed a combined MCMC analysis incorporat- 
ing the DECPHOT 8 |i.m hght curves, which were detrended 
with Equation |3] The raw and detrended data are shown 
with the best-fitting trend and occultation models in Fig- 
ure |4] We found consistent 8 \im occultation depths when 
incorporating only one dataset or both datasets in our anal- 
ysis (Table O). The residuals of the decphot light curves 
exhibit a slightly smaller scatter than the aperture photom- 
etry light curves do (Figure ID Table [3)1 . These decphot 
depths and associated uncertainties are in close agreement 
with those derived using the hght curves obtained from sim- 
ple aperture photometry (Table [3}. This is also the case for 
ecoso;, esinoj and the time of mid-occultation (Table IS}. 
Thus our method of obtaining 8 [im hght curves by simple 
aperture photometry and detrending them with detector po- 
sition is verified, and it is these hght curves that we use in 
the simultaneous analysis from which we calculated our sys- 
tem parameter values. 
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3.4 Partitioning of data 

In our simultaneous MCMC analysis we partitioned the 
WASP photometry according to observation season and 
camera into five datasets, so that each dataset could thus be 
normalised independently, as was done in AlO. As in TIO, 
we partitioned the RV data into four datasets: CORALIE 
data sampling the full orbit (rvl)\ HARPS data sampling 
the full orbit (rv2); a spectroscopic transit and a week of 
adjoining data as measured by CORALIE (ru3); a spectro- 
scopic transit and two RVs from the following day as mea- 
sured by HARPS (rw4). Both an instrumental offset and a 
specific stellar activity level have the potential to affect the 
measured RV of a star. The spectroscopic transits comprise 
a large number of RVs taken in quick succession, whereas 
the data sampling the full orbit were taken over a long time 
span and are thus expected to sample a range of stellar activ- 
ity level that should average to a mean value of zero (TIO). 
Thus, by partioning the RV data, we allow each dataset to 
have its own centre-of-mass velocity 7, thus avoiding the 
risk of obtaining spurious values for the planet's mass and 
orbital eccentricity. 



3.5 Photometric and RV noise 

We scaled the photometric error bars so as to obtain a re- 
duced of unity, applying one scale factor per dataset. The 
aim was to properly weight each dataset in the simultane- 
ous MCMC analysis and to obtain realistic uncertainties. For 
the five sets of WASP photometry the scale factors were in 
the range 0.87-0.96. The error bars of the Euler photometry 
were multiplied by 1.33. The scale factors for the occulta- 
tion photometry were in the range 1.04-1.09. Importantly, 
the error bars of the occultation photometry were not scaled 
when deciding which trend models or aperture radii to use. 

We assessed the presence of correlated noise in the 
Spttzer and Euler data by plotting the RMS of their binned 
residuals (Figure [7|. Though there is no correlated noise ev- 
ident in the Spitzer data, it is present at a small level in 
the Euler data over time-scales of 8-80 minutes. Due to the 
similarity with the time-scales of the fitted features in the 
transit (ingress takes 36 minutes, as does egress, and the 
transit duration is 264 minutes), the values of some fitted 
parameters may be affected to a small degree. 

For the same reasons as with the photometry, we added 
a jitter term in quadrature to the formal radial velocity er- 
rors, as might arise from stellar activity. We used an initial 
MCMC run to determine the level of jitter required for each 
dataset to obtain a reduced of unity. We found that the 
HARPS orbital data (r«2) required a jitter of 3 ms~^ and 
the HARPS spectroscopic transit data (rvi) required a jit- 
ter of 20 ms~^. It was not necessary to add any jitter to 
either of the two CORALIE datasets. 



3.6 Time systems and light travel time 

The Euler photometry and the CORALIE RVs are in the 
BJD (UTC) time system. The WASP and Spitzer photom- 
etry are in the HJD (UTC) time system. The difference be- 
tween BJD and HJD is less than 4 s and so is negligible for 
our purposes. Although leap second adjustments are made 
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Figure 7. RMS of the binned residuals for the new Spitzer occul- 
tation photometry (upper four panels, with the datasets presented 
in the same order as in Figurc[T]l and the existing Euler photome- 
try (lower panel). The solid black lines, which are the RMS of the 
unbinned data scaled by the square root of the number of points 
in each bin, show the white-noise expectation. The ranges of bin 
widths (1-250 minutes for Spitzer and 2.5-180 minutes for Euler) 
are appropriate for the datasets' cadences and durations. 



to the UTC system to keep it close to mean solar time, mean- 
ing one should really use Terrestrial Time, our observations 
span a short baseline (2006-2008), during which there were 
no leap second adjustments. 

The occultation of WASP-17b occurs farther away from 
us than its transit does, so we made a first order correction 
for the light travel time. We calculated the light travel time 
between the beginning of occultation ingress and the begin- 
ning of transit ingress to be 50.4 s. We subtracted this from 
the mid-exposure times of the Spitzer occultation photome- 
try. As we measure the time of mid-occultation to a precision 
of ±150 s, the impact of this correction was small. 



4 RESULTS 

TableOshows the median values and the 1-a uncertainties of 
the fitted proposal parameters and derived parameters from 
our final MCMC analysis. Figure [1] shows the best-fitting 
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trend and occultation models together with the raw and de- 
trended Spitzer data. Table [2] gives the best-fitting values 
for the parameters of the trend models (Equation [2]), to- 
gether with their 1 a uncertainties. Figure |8] displays all the 
photometry and RVs used in the MCMC analysis, with the 
best-fitting eclipse and radial velocity models superimposed. 

From this we see that WASP-17b is a very bloated 
planet (_Rpi — 2.0 -Rjup) in a slightly eccentric, 3.7 day, ret- 
rograde orbit around an F6V star. By constraining the ec- 
centricity of WASP-17b's orbit to low values we have shown 
that the circular solution presented in AlO (in which a total 
of three solutions were presented) is closest to reality. 



4.1 Orbital eccentricity 

We have shown the orbit of WASP-17b to be non- 
circular: ecoscj is non-zero at the 4.8-a level (ecoso; — 
0.00352j;o 00073; Figure[ni), and the best-fitting solution sug- 
gests that WASP-17b is occulted by its host star 12.0 ± 
2.5 min later than if it were in a circular orbit. Our mea- 
surement of e cos uj rules out large values of e for all orbital 
orientations other than those with |cj| ~ 90, and the lim- 
its we place on esinoj prohibits large values of e for those 
orientations with \uj\ ^ 90 (Figure [TO)) . From the MCMC 
analysis, the 1-cr (68.3 per cent) lower and upper limits on 
e are, respectively, 0.010 and 0.043 and the 3-cr (99.7 per 
cent) lower and upper limits on e are, respectively, 0.0019 
and 0.0701. We can set a more stringent 3 a lower limit on 
e by assuming esinoj ~ (and so ~ 90), in which case 
it would be equal to that of the 3-a lower limit on ecostj: 
0.0012. 

Almost all values of a; are permitted by the current data, 
with only ~ 90 being ruled out by the limits placed on 
e sin ij (Figure llOf) . Large values of e are consistent with the 
data only if ~ 90, otherwise any orientation of the orbital 
major axis is permitted providing that e is small. We can 
thus use our measurement of e cos u to infer a probable value 
of e. For random orientations of the major axis, the expected 
value of cosoj is E(cosaj) — 2/n. Thus, the expected value of 
e is E(e) = e coscj/ E(cos cj) = 0.0055. 

We explored the effect of each occultation photometry 
dataset in turn on the orbital eccentricity, and of all four 
datasets combined. We did so by performing MCMC runs 
that incoporated either all, none or just one of the Spitzer 
datasets (Table (6] Figure (9)1 . This demonstrates how valu- 
able the Spitzer occultation photometry is in determining 
orbital eccentricity, as its inclusion in our combined anal- 
ysis caused the size of the 68.3 per cent confidence inter- 
val for ecostj to decrease by a factor 30.2, and the inter- 
val for esintj to decrease in size by a factor 2.4. In addi- 
tion to the RV data, it is the orbital phase of the occulta- 
tion that constrains ecosoj and it is the occultation dura- 
tion, relative to the transi t duration, that constrains esino) 
l|Charbonneau et al.ll2005l ). When including any one of the 
four occultation datasets, the best-fitting values of ecosoj 
and e sin lj obtained are consistent with the values obtained 
when including all four datasets. Thus no individual dataset 
is biasing our best-fitting solution. 
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Figure 9. Top panel: A comparison of the posterior probabil- 
ity ciistributions of e cos uj and e sin uj from our combined MCMC 
analysis when including (red dots) and excluding (cyan dots) the 
occultation photometry. The extent of the error bars show the 
1-cr confidence limits and their intersections show the median val- 
ues. Middle panel: Normalised histogram of the ecosw posterior 
probability distribution from our combined MCMC analysis in- 
coroporating the Spitzer photometry. The point with error bars, 
arbitrarily placed at probability = 0.01, depicts the best-fitting 
value and its l-cr error bars. Bottom panel: The same plot as the 
middle panel, but for esinoj. 
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Table 5. System parameters of WASP- 17 



Parameter 


Symbol 


Value 


Unit 


Orbital period 


P 


3.7354380 ± 0.0000068 


d 


Epoch of mid-transit (HJD, UTC) 


Tc 


2454577.85806 ± 0.00027 


d 


Transit duration 




0.1830 ±0.0017 


d 


Duration of transit ingress Rd duration of transit egress 


Ti2 ~ T34 


0.0247 ±0.0017 


d 


Planet / star area ratio 




0.01696 ± 0.00026 




Impact parameter 


b 


n 401+0 059 

U.4Ui_Q Qyy 




Orbital inclination 


i 


86.83l°;i| 





O liciidl 1 diiidl 1 ciitiJV V^iUClLV 


^ 1 




ms-1 


Semi-major axis 


CL 


n ^^^^ -t- n nnn'^4 


ATT 


Centre-of-mass velocity 


O'rvl 


-49 513.67 ± 0.56 




Offset between RV dataset rv2 and rvl 


7rv2-rvl 


22.07 ±0.68 




Offset between RV dataset rv3 and rvl 


O'rvS — rvl 


13.5 ±2.2 


ms"-'- 


Offset between RV dataset rv4 and rvl 


0'rv4 — rvl 


25.6 ±2.8 




Orbital eccentricity 


e 


0.028±0;Ooil 








0.0019 < e < 0.0701 (3 a) 




Expectation value of orbital eccentricity 


<e> 


0.0055 




Argument of periastron 




-82.611" 







e cos Lo 


f\ n,^or--n-l-0 00076 

0.0035210 00076 








-0 027+0-019 

"■'-'^'-0.015 




Phase of mid-occultaion, having accounted for light travel time 


'^^rnid — o cc . 


0.50224 ± 0.00050 




Occultation duration 


T58 


0.1746+°Ooil 


d 


Duration of occultation ingress fa duration of occultation egress 




0.0232 zb 0.0016 


d 


T-? Ci \ 'i" 1 "\ m T^ltlTlO'l" "Hll'V f^'i" A- ^ 11 Tn 

J.\,t;iajLi VU JJlttliUL aLOji iiUJl. cLLr '-t.xj [^ili 


^^4.5 f^m 


00990 -1- 0001 




Relative planet-star flux at 8 ^.m 


AFg i^m 


0.00237 ± 0.00039 




Planet brightness temperature! at 4.5 \ivci 


^b,4.5|a.m 


1881 ± 50 


K 


Planet brightness temperature^ at 8 |J.m 


Tb,8|j.m 


1580 ± 150 


K 


Sky-projected stellar rotation velocity 


V sin / 




km s~ 


Sky-projected angle between stellar spin and planetary orbit axes 


A 


-148.71^;^ 





Star mass 


M* 


1.306 ± 0.026 




Star radius 


R, 


1.572 ±0.056 


Rq 


Star density 


p. 


0.336 ± 0.030 


Pe 


Star surface gravity 


logg. 


4.161 ±0.026 


(cgs) 


Star effective temperature 


Teff 


6650 ± 80 


K 


Star metallicity 


[Fe/H] 


-0.19 ±0.09 




Planet mass 


Mpi 


0.486 ± 0.032 




Planet radius 


^pi 


1.991 ±0.081 


Rjup 


Planet density 


Ppi 


0.0616 ±0.0080 


Pjup 


Planet surface gravity 


loggp 


2.448 ± 0.042 


(cgs) 


Planet equilibrium temperature^ (full redistribution) 


'^P,A=0,f=l 


1771 ± 35 


K 


Planet equilibrium temperature^ (day side redistribution) 


Tp,A=0,f=2 


2106 ±41 


K 



t We modelled both star and planet as black bodies and took account of only the occultation depth uncertainty, which dominates, 
t - fit 

f Tp A=o,f = /^"^cffy "2^ where / is the redistribution factor, with / = 1 for full redistribution and / = 2 for day side redistribution. 
We assumed the planet albedo to be zero, A = 0. 



5 DISCUSSION 



5.1 Planet radius 



With a radius of 2.0 -Rjup, WASP-17b is the largest known 
planet by a margin of 0.2 -Rjup, and is over 0.7 -Rjup larger 
than predicte d by standard cooling theory of i rradiated gas 

gia nt planets ( Fortnev. Marlev fc Earned [20071). 

[ ibgui fc Burrows! (|2009l ) and llbgui. Spiegel fc Burrows! 
(|201ir used a coupled radius-orbit evolutionary model to 



show that planet radii can be inflated to 2 -Rjup and beyond 
during a transient phase of heating caused by tidal circular- 
isation of a short (a ~ 0.1), highly eccentric (e « 0.8) orbit. 
Though, as was noted in both studies, planets can persist in 
an inflated state for Gyr beyond the circularisation of their 
orbit and the cessation of tidal heating, they do cool and 
contract significantly prior to full circularisation. In each 
study the orbits are still significantly non-circular (e > 0.1) 
when the planets are largest. Thus, under the transient heat- 
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Figure 8. The results of our combined analysis, which combines the new Spitzer occultation photometry with existing photometry 
and radial velocity measurements. The models generated from the best-fitting parameter values of Table [5] are overplotted. Top-left: 
Occultations at 4.5 |im and, offset in relative flux by —0.01, 8 [xm. The two occultations per channel from Figure [l] were binned 
(A</> = 0.002, ~11 min) together. Top-right: Transit light curve taken with Euler in the /c-band (data from AlO). Middle: Photometric 
orbit and transit illustrated by WASP-South data (data from AlO). Bottom: Spectroscopic orbit and transit illustrated by CORALIE 
and HARPS data (data from AlO and TIO). The measured systemic velocities of each dataset (Table [Sj have been subtracted. 



Table 6. Effect of occultation light curves on best-fitting orbital eccentricity 



Included occultation photometry 


e 




ecoso; 


esinw 


Tocc T^occ, circular (min)^ 


4.5 ^m, 2009 Apr 24 
4.5 nm, 2009 May 1 
8 tim, 2009 Apr 24 
8 n-m, 2009 May 1 
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All (decon. phot.) 
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0055+° °°^^ 
U.UUOO_o 0024 

U.U10_o 0J2 
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0.0381°:°^^ 
0.028tHi^ 
0.022tni^ 


1 0+82 
-^■^-56 
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00302+° °°l°^ 
U.UUOUZ_o 00098 

U.UUZi_0 0067 
00662+° °°°°^ 

U.UUDOZ_o 00111 

011+°-°27 
U.Uii_o 018 

00352+° °°°'^'' 
U.UUOOZ_o 00073 

n0335+° °0°'^3 
U.UUcWO_o 00075 


-0 nt;i+0 020 
U.UOi_o 017 

-0 001+°"°* 
U.UUl_o 007 

-0 005+° °" 
U.UUO_o 062 

-0 049+° °2* 
U.U'i»_o 020 

r, (,,0-1-0.062 

U.UiJ_o 021 

-0 027+° °^^ 
'J-'-'^' -0.015 

-0 022+° °" 
U.UZZ_o 016 


12.7 ± 2.9 

10.3^3.4 

7 q-l-13.2 

'•■^-23.2 
22.7tli 

-61 
12.0 ±2.5 
11.5 ±2.6 



^ Tocc is the time of mid-occultation derived from a simultaneous MCMC analysis. 



ing scenario, the very largest planets are expected to have a 
non-zero eccentricity. Though we do measure a non-zero ec- 
centricity for WASP-17b, it is small, and the stringent upper 
limit that we place on e is inconsistent with current models 
of one transient phase of tidal heating. 

Other than transient heating, ongoing tidal heat- 
ing may occur if the orbit of a planet were kept non- 
circular by the continuing inte raction with a third body 
jlbgui. Burrows fc SpiegefcoiOl ). However, the stringent up- 
per limit we place on e makes this unlikely as the sole 



cause of the inflation of WASP-17b, as it would necessi- 
tate a lower planetary tidal dissipation factor than theo- 
retical models or empirical deter minations generally suggest 
(|lbgui. Burrows fc SpiegelbOlOl ). 

If the atmospheric opacity of WASP-17b were enhanced 
then its internal heat woul d be lost at a lower r ate and con- 
traction would be slowed (|Burrows et al.|[2007h . The atmo- 
spheric opacities of WASP-17b may be enchaned if, for ex- 
ample, the strong optical and UV irradiation of the planet by 
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Figure 10. The range of e and u) permitted by the available 
data. The black dot with error bars shows the best-fitting values 
from our combined MCMC analysis. The grey dots are the values 
in accepted MCMC steps. The solid blue and red lines show the 
values of e and oj that would be indicated by, respectively, the 
best-fitting values of e cos ui and e sin u) on their own. The dashed 
and the dotted lines bound the parameter space permitted by the 
l-cr and the 2-cr limits, respectively, on those parameters, with the 
same colour scheme applying. Note that, as the 2-a upper limit 
on esint<; is positive, almost all values of u) are consistent with 
the data at the 2-cr level (providing e < 0.01). 

its host star produces thick hazes, absorbing clouds and non- 
equilibrium chemical species (e.g. tholins or polyacetylenes) . 

The bloated planets are all very strongly irradiated by 
their host stars, and a small fraction of stellar insolation en- 
ergy would be sufBcient to accou n t for the observed degrees 
of bloating. iGuillot fc ShowmanI l|2002l) suggested that the 
kinetic energy of strong winds, induced in the atmosphere 
of a short-period planet by the large day-night temperature 
contrasts that result from tidal locking, may be transported 
downward and deposited as thermal energy in the deep in- 
terior. However, a mechanism to convert th e kinetic energy 
into t herm al energy wo uld still be required. iLi fc GoodmanI 
l|201(]| ) and lYoudin fc M itchell (2010) found that turbulence 
is efficient at dissipat ing kinetic energy. Magnetic dra g on 
weakly ionized winds (jPerna. Menou fc Rauscheill2010l ) and 
Ohmic heating (jBatygin fc StevensoDll201Gl ) are alternative 
mechanisms. 

5.2 Planetary atmosphere 

iFortnev et al.l l|2008l ) hypothesise that the presence of high 
opacity TiO and VO gases in the atmospheres of highly irra- 
diated planets (those experiencing an incident flux of > lO'' 
erg s~^ cm~'^) cause them to have temperature inversions. 
Thus, with an incident flux of 2.2 ± 0.2 x 10® erg s~^ cm~^, 
WASP- 17b is expected to have an atmospheric temperature 

inversion unde r this hypothesis. 

However, ISpiegel. Silverio fc Burrow^ (|2009l ) suggest 
that, for a planet with the insolation level of WASP-17b, 
it is unlikely that a temperature inversion could be caused 
by the presence of TiO and VO in the upper atmosphere. 
They find that a cold trap exists between the hot convection 
zone and the hot upper atmosphere on the irradiated day 



side, in which titanium is likely to form condensates that 
settle more strongly than does gasesous TiO. Therefore, un- 
less there is extremely vigorous macroscopic mixing and the 
condensed Ti is lofted back in to the upper atmosphere then 
it is unlikely that TiO can explain the observed temperature 
inversions. Not only does VO have the same 'cold trap' is- 
sue, but it also has a lower opacity than TiO and is an order 
of magnitude less a bundant. 

[Knutson. Howard fc Isaacson! (|20ld ) suggest that plan- 
ets orbiting chromospherically active stars do not have tem- 
perature inversions, and planets orbiting quieter stars do 
have inversions. They suggest that the high UV flux that 
planets orbiting active stars are likely to experience de- 
stroys the compo unds responsible for the obser v ed tem per- 
ature inversions. iKnutson. Howard fc Isaacson! l|2010l ) find 
the two classes to be delineated by a host-star activity level 
of log(_RHK) ~ —4.9. Though they caution that the calibra- 
tion for log(_Rjjj;) is uncertain for stars as hot as WASP-17, 
they measure the star to be quiet: log(_RHK) = —5.3. WASP- 
17b is therefore expected to have an atmospheric tempera- 
ture inversion under this hypothesis as well. 

In Figure 1111 the measured 4.5 and 8 \x.m planet-star 
fiux-density ratios are compared to two model atmospher e 
spectra of the planet (|Barman. Hauschildt fc AUardllioOsI ). 
with parameters taken from TableO A black-body (Tp_A=o= 
1600 K) is a poor fit to the data and is thus ruled out. In one 
model atmosphere TiO produces a temperature inversion 
across the photospheric depths. In the other model, there 
is no atmospheric TiO. The two models have near-identical 
4.5 and 8 \ira absolute fluxes, and so we can not currently 
discriminate between the two. A precise measurement at 3.6 
\im may distinguish between the two cases and thus reveal 
whether WASP-17b has an atmospheric temperature inver- 
sion. 

By modelling the planet and star as black bodies, we 
used the measured planet-star flux-density ratios to calcu- 
late 4.5 and 8 \x.m brightness temperatures of 1881 ± 50 K 
and 1580 ± 150 K, respectively. We calculate an equilibrium 
temperature rp,A=o,f=i = 1771 ± 35 K by modelling the 
planet as a black body with efficient redistribution of en- 
ergy from its day side to its night side. The closeness of the 
brightness temperatures to this equilibrium temperature is 
consistent with the planet having a low albedo and efficient 
heat redistribution. 



5.3 Misaligned orbit 

WASP-17b is in a retrograde orbit. For planet-planet or star- 
planet scattering to have caused the misalignment between 
the orbit of WASP-17b and the spin axis of its host star, 
an additional body must have been present. We looked for 
evidence of a long-term drift 7 in the radial velocity mea- 
surements, which span 716 days, as may be caused by the 
presence of a long-period companion. From a straight-line 
fit to the residuals of the radial velocities about the best- 
fitting model, we get 7 = —6 ± 5 ms~^ yr"'^- Hence, there 
is currently no evidence for a third body in the system, but 
this does not preclu de planet-planet scattering as th e cause 
of the misalignment. iNagasawa. Ida fc Besshd (|2008l ) found, 
whilst showing that a combination of planet-planet scatter- 
ing and the Kozai mechanism can put planets into short. 
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Figure 11. Comparison of planet-star flux density measurements 
with two model planet atmospheres and with a black-body. The 
model atmosphere with TiO exhibits a temperature inversion that 
extends down to photospheric depths, whilst the model without 
TiO does not. Inset: Temperature-pressure profiles for the two 
model atmospheres. 
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Figure 12. Modified H-R diagram. The isochrones (Z = 0.012 ^ 
[Fe/H] = -0.19) for the ag es 0.5, 1, 2, 2.5, 3, 4, and 5 Gyr 
are from iMarigo et al.l l|2008l ) and the evolutionar y mass tracks 
Z = 0.012 ^ [Fe/H] = -0.019; Y = 0.30) are from lBertelli et alj 
20081 ). To obtain the mass tracks, we performed a simple linear 
interpolation of their Z = 0.0008 and Z = 0.017 tracks. 



lithium abundance (Ali < 1.3) in AlO is consistent with 
this. Thus, lithium is not an effective indicator of age for 
WASP-17. 



retrograde orbits, that the outer planets can end up at large 
orbital distances, making them difficult to detect, or they 
can be ejected from the system. 

5.4 System age 

We in terpolated the stellar evolution tracks of lMarigo et al.l 
using p, from Table [5] and the values of T^a and 
[Fe/H] from TIO (Figure [T2)) . This suggests an age of 2.65 ± 
0.25 Gyr and a mass of 1.20 ± 0.05 Mq for WASP-17. 

Assuming the stellar-spin axis to be in the sky plane, the 
measured «sin/ of WASP-17 and its derived stellar radius 
(Table [S]) indicate an upper limit to the rotational period 
of Prot = 7.91 ± 0.75 d. C ombin i ng th is with the B ~ V 
colour of an F6V star from iGravl (|2008h . and the relation- 
ship of iBarned l|2007l ). we estimate an upper limit on the 
gyrochronological age of 1.9 ± 0.5 Gyr. We found no evi- 
dence for rotational modulation in the WASP light curves. 

We calculated a tidal circularisation time-scale of Tcirc = 
5 Myr for WASP-17b by using the best-fitting values of the 
planetary (Qp = lO'^'^) and stellar (Q^ = 10^-^) tidal dis- 
sipation factors of Ijackson. Greenberg fc Barnes! (fioOSa) in 
their Equation 1. As the values of the tidal dissipation fac- 
tor are highly uncertain (Qp = lO'^-lO*, Q, = 10^-10*, e.g. 
Ilbgui. Spiegel fc Burrowil201ll ). a range of Tdrc = 2-1700 
Myr is possible. 

With Teft = 6650 ± 80 K (TIO), WASP-17 is in the 
'Lithium gap' (or 'dip'), which is the range of T^fi — 6600 ± 
150 K in which stars are depleted in lithium b y a factor of 30 
ormore than in hotter and cooler stars fsee iBalachandranI 
1 19951 . and references therein). The upper limit placed on the 



6 CONCLUSIONS 
6.1 Science 

With a radius of 2.0 -Rjup, WASP-17b is larger than any 
other known planet by 0.2 -Rjup and it is 0.7 i?jup larger 
than predicted by standard cooling theory of irradiated gas 
giant planets. The extent of the planet's inflation is difficult 
to explain with current models. 

Our Spitzer occultation photometry gives much tighter 
constraints on orbital eccentricity than existing radial veloc- 
ity data alone, thus permitting an accurate determination of 
the stellar and planetary radii. We have shown that WASP- 
17b is in a slightly eccentric orbit, with 0.0017 < e < 0.0701. 
The stringent upper limit we have placed on eccentricity sug- 
gests that a transient phase of tidal heating alone could not 
have infiated the planet to its measured radius. Nor could 
ongoing tidal heating involving a third body, unless the plan- 
etary tidal quality factor is smaller than the best theoretical 
and empirical determinations. 

We find no evidence in the radial velocity measurements 
for a third body in the system, the presence of which would 
be necessary to excite the eccentricity of WASP-17b for tidal 
heating to be ongoing, and may have been necessary to mis- 
align the planet's orbital axis with the spin axis of the star. 

Our 4.5 and 8 (xm planet-star flux-density ratios do not 
probe the existence of the expected atmospheric tempera- 
ture inversion, but a measurement at 3.6 [im may do so. 
Though the ratios are inconsistent with a black-body atmo- 
sphere, they are consistent with a low-albedo planet that 
efficiently redistributes heat from its day side to its night 
side. 
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6.2 Spitzer data 

To determine correctly the photometric uncertainties and 
the optimal aperture radii to use for Spitzer data, account 
must be taken of the counts removed during sky-dark sub- 
traction. 

When the background is bright relative to the target 
at 8 \im, the measured occultation depth can depend sen- 
sitively on the choice of aperture radius. In these circum- 
stances detrending with detector position vastly reduces the 
dependency. An alternative is to perform deconvolution pho- 
tometry. 

In addition to the known hour-long oscillations of 
Spttzer's pointing about the nominal position, there is also 
a high-frequency jitter, with periods of 1-2 minutes. So, 
when accounting for the inhomogeneous detector response 
(or 'pixel phase' effect), one should detrend target flux with 
the unsmoothed target detector positions. 
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